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ABSTRACT 

The production rate of ionizing photons in young (< 8 Myr), unresolved stellar clusters in the 
nearby irregular galaxy NGC 4214 is probed using multi-wavelength Hubble Space Telescope WFC3 
data. We normalize the ionizing photon rate by the cluster mass to investigate the upper end of the 
stellar initial mass function (IMF). We have found that within the uncertainties the upper end of 
the stellar IMF appears to be universal in this galaxy, and that deviations from a universal IMF can 
be attributed to stochastic sampling of stars in clusters with masses ^ 10"^ Mq. Furthermore, we 
have found that there does not seem to be a dependence of the maximum stellar mass on the cluster 
mass. We have also found that for massive clusters, feedback may cause an underrepresentation in 
Ha luminosities, which needs to be taken into account when conducting this type of analysis. 

Subject headings: galaxies: individual (NGC4214) - galaxies: star clusters: general - galaxies: star 
formation - stars: luminosity function, mass function - stars: massive 



1. INTRODUCTION 

The stellar initial mass function (IMF), the distribu- 
tion of stellar masses in newly formed stellar popula- 
tions is essential for understanding the evolution and 
star formation histories of galaxies. Whether it is uni- 
versal or dependent on environment has been a highly 
contested issue over the past few years. While IMF 
measurements in high density environments like the 
Milky Way and Magellanic Clouds have indicated an 
invar iant IMF (jBastian et al.l 120101: iMassevI 120031: lOevI 
1201 It ) , other claims of a non-universal IMF have been 
made (jvan Dokkum fc ConrovllMlt iWilkins et aT1[200l 
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Fardal et al.1 120071 : iKroupa et all 120111 : iCappellari et al.l 
2012f ). In particular, star-forming dwarf galaxies may 



show a deficiency in the ionizing photon rate per unit 
UV or optical luminosity (Hovcrsten fc Glazcbroo k 200g 



Lee et al. 2009; Meurer et all 120091: iBoselli et allfeOOS 
Gunawardhana et al.. ■2011 ). However. iFumagalli et all 



( 201 It ) and lWeisz et al.l ()2012f ) have shown that stochas- 
ticity in populating the IMF or bursts of star formation 
can explain the observed variations in IjHa/^FUV so that 
an unusual IMF is not required. Clearly this is an unre- 
solved issue that is essential to understanding the funda- 
mental evolution of galaxies. 

Traditional methods for IMF measurements in the 
Milky Way, Large Magellanic Cloud (LMC), and Small 
Magellanic Cloud (SMC) are to count individual stars 
in clusters < 3-5 Myr old tha t still retai n their most 
massive stars (lAnderson et al.l 120091 : ISabbi et al.l l2008i : 
iSirianni et al.ll2000l for example) . Generally only 20% of 
stellar cluster s will survive early mass loss to live longer 
than 10 Myr (jLada fc Ladall2003l . "infant mortality"), so 
catching them very early is essential for observing the 
full stellar population. Even nearby, significant crowd- 
ing can cause these star counts to be incomplete, with 
high mass stars and low mass stars suffering from differ- 
ent selection biases: low mass stars are generally harder 
to count, due to the inability to easily detect smaller, 
fainter stars and due to dynamical ejection of the low 
mass stars, while massive stars may suffer from confu- 
sion due to their sinking towards the cente r of the cluster 
()Ascenso et al.l2009HMaiz Apellani"gi2008D . This method 
becomes progressively less effective at distances outside 
of the Magellanic Clouds (~ 50 kpc), even with the res- 
olution of the Hubble S pace Telescope (HS T). 

As demonstrated in Calzetti et al.l ()2010l ) , it is possible 
to constrain the upper end of the IMF in external galax- 
ies without the use of individual star counts. Instead, 
a nearly coeval population that still contains the most 
massive members capable of producing ionizing photons 
can be constructed from the sum of unresolved young 
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Fig. 1.— Color composite WFC3 image of NGC 4214 courtesy of Zoltan Levay (STScI-2011-14), R. O'Connell (GO 11360), and the 
WFC3 SOC. The footprint of the WFC3/IR images is shown in white. Blue circles are clusters with undetected Ha emission, yellow circles 
are clusters with Ha emission. There are a total of 52 compact clusters. 



clusters in the galaxy. Therefore, measuring Q(H"), the 
hydrogen ionizing photon rate, in these young clusters 
is equivalent to measuring the number of massive stars. 
This method relies on normalizing the ionizing photon 
rate to the age-independent cluster mass, which is anex- 
tension of the method described in Cor belli et al.l (j2009D . 
The treatment presented here can eliminate the need for 
age-dependent bolometric luminosities, but does require 
that the cluster ages are ac curately determined. The pi- 
lot study done on M51a bv lCalzetti et al.l (|2010l ). found 
that there was no obvious dependence of the upper mass 
end of the IMF on the mass of the star cluster down 
to - 10^ Mq, but a more extensive analysis including 
additional populations than M51a is needed for conclu- 
sive evidence. Specifically galaxies with star formation 
rates (SFR) below the threshold for which IMF variances 
have been suggested (< 0.1 M© yr^^) need to be inves- 
tigated. This paper aims to extend this study using a 
nearby galaxy with a lower SFR more similar to those 



dwarf galaxies which may be exhibiting a deficiency of 
ionizing photons. 

NGC4214 is an irregular, LMC-type, star bu rsting 
galaxy located ~3 Mpc away (jPopita et al.ll2010l ) with 
an Ha and UV SFR of 16 Mp yr-^ and 0.22 M© 
yr"^ respectivel y (|Lee et al.ll2"0lll:lKennicutt et al.ll2008l: 
ILee et ail 120091) as weU as a sub-solar (Z^ 0.25 Zg) 
metallicitv ()Kobulnickv fc Skillmanlll996[) . The star for- 
mation history (SFH) of the central region of NGC4214 
shows a strong increase in SFR starting 100 Myr ago 
with a prominent peak at recent times (< 10 Myr); de- 
spite this, less than 1% of the mass of the galaxy is due t o 
the current star formation event ([Williams et al.l [20 111 ). 
Due to its proximity and recent star forming activity, 
NGC 4214 is an ideal test-bed for the upper end of the 
IMF (uIMF). In Section 2 of this paper we will discuss 
the observations and cluster selection criteria, in Section 
3 we will present the models and age and mass determi- 
nations, and in Section 4 we discuss the results. 
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2. OBSERVATIONS AND CLUSTER SELECTION 

Observations were taken with HST WFC3/UVIS 
and WFC3/IR as part of GO 11360 (PI: O'Connell). 
The observations on which we concentrate here in- 
clude F225W (1665s), F336W (1683s), F438W (1530s), 
F547M (1820s), F657N (1592s), F814W (1339s), F HOW 
(1198s), and F128N (1198s), shown in FigurelH For ease 
of discussion, we will refer to these as NUV, U, B, V, Ha, 
I, J, and P/3 respectively. Each flat-fielded image was 
co-added, cosmic rays were removed, and corrections for 
distortion were made using the task MULTIDRIZZLE 
into a final pixel scale of 0".0396 pixel"^. This corre- 
sponds to 0.5 8 pc pixel"^ at a distance of 3 Mpc. See 
IDopita et al.l ([2010) for a full explanation of the reduc- 
tion procedure, including the creation of continuum sub- 
tracted Ha images. A similar procedure was adopted to 
create a continuum subtracted P/3 image, which covers 
~75% of the UVIS field of view. The footprint of the IR 
field of view is also shown as a white outline in Figure [TJ 

Cluster candidates were identifie d using a technique 
similar to that described in Chanda r et al.l (j2011l 120101) 
for M51 and M83. Aperture photometry was performed 
on the wide and medium band images using the IRAF 
task PHOT with an aperture of 3 pixels, and a back- 
ground annulus between 10 to 13 pixels. The aperture 
corrections were addressed similarly as was done on clus- 
ters in M83 from [Chandar et al.l ()2010D . In their study 
they chose two different methods for aperture correction, 
both of which rely on the concentration index (C, the 
difference in magnitudes between 3 pixel and 0.5 pixel 
radius). Method one used a single value for the aper- 
ture corrections of point sources (C < 2.3) and a slightly 
larger value for extended sources (C > 2.3). For this 
sample we chose to use their second approach, which is to 
use an aperture correction equation for extended objects 
with 2.3 <C< 3.4. Photometric conversion from counts 
to erg cm~^ s~^ were accomplished using the filter 
dependent PHOTFLAM values provided by the STScI 
website. Galactic foreground extinction of E{B — V) = 
0.02 was correcte d using the Milky Way extinction curve 
from lFitzpatrickl (|l999). 

Due to the more extended nature of HII regions 
surrounding the stellar clusters, aperture sizes that 
scaled with the cluster mass according to the expected 
Stromgren radius were used to measure the hydrogen re- 
combination lines on the continuu m-subtracted Ha+ [N 
II] and P/3 images. As was done in lCalzetti et al.l ()2010[ ). 
a radius of about 0.35 ^stromgren was selected due to the 
crowding of the clusters, which corresponds to a range 
between 6-10 pixels. This radius is sufficiently larger 
than the PSFs for both the Ha and Pf3 images, so there 
are no concerns of PSF variations. The local background 
was subtracted using a 3 pixel wide annulus centered 
on the cluster outside of the aperture radius in order to 
avoid contamination from other diffuse emission. Aper- 
ture corrections were calculated from a few, very iso- 
lated sources, and were applied to the other regions. We 
found that an additional correction of 1.30 times the flux 
was needed to go from a 20 pixel aperture to an "infi- 
nite" aperture. Contamination from [N II] was removed 
using the average galact ic [N II]/Hq; ratio of 0.11 from 
[Kobulnicky & Skillman (|1996D . The corrections for the 
ionized gas extinctions were measured region by region 



using the WFC3/IR P/3 image for those regions where 
the IR a nd UVIS overla p (see Figure [B using the formu- 
lation in 'Cal zetti et al.l (|2000f) and were applied to the 
Ha luminosities. 

The cluster catalog was populated by the combina- 
tion of two methods. The first, and most robust, 
was using the a utoma ted method discussed in depth in 
I Chandar et al.l (|2010[ ) which accounts for roughly 80% 
of the cluster catalog. To make sure all clusters are ac- 
counted for, we also use a manual selection procedure 
from a careful examination of the WFC3 images. This 
ensures we identify clusters in crowded regions or clus- 
ters near a bright star which may have been missed in 
the automated process. In total we have identified ^ 400 
cluster candidates. 

3. ANALYSIS 

In this section we will present an in-depth discussion 
of the analysis procedure used in this paper, and a de- 
tailed description of the filter convolution is included 
in the Appendix. For the mass and age determina- 
tion of the clusters we have used broad-band photom- 
etry (without the Ha included) to determine ages and 
masses of the clusters using both a canonical and trun- 
cated IMF from stochastic and deterministic stellar mod- 
els. The best fit model-derived masses of clusters with 
ages less than 8 Myr are then binned into three dis- 
tinct mass bins. Within these mass bins the masses 
are summed and the Ha luminosity are summed to de- 
termine the Ljjq/Mc/ ratios for various cluster masses 

"iff X]' ^' These ratios are then compared to 
predicted models with two different assumptions about 
the I MF, the canonical one and the variable upper mass 
limit (|Kroupa fc Weidneiil2003l . see Discussion), in order 
to constrain the IMF of NGC 4214. 




Age [Myr] 

Fig. 2. — Comparison of V-band luminosities of 1 X 10^ Mq 
SLUG models from 1-20 Myr with SB99 (red line) with Av' = 
0. The magnitude of scatter from the stochastic models will also 
create a scatter of cluster mass. Note that the arithmetic mean of 
the SLUG models match the SB99 models, but that the geometric 
mean is slightly lower. 

3.1. The Models 

In order to accurately determine ages and masses of 
the clusters, the spectral energy distributions (SEDs) 
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created from the photometry of the NUV, U, B, V, I ob- 
servations were compared to stellar synthesis models. 
While the norm has been to use deter ministic model s 
such as STARBURST99 (SB99) (Leithe rer et al.l[T999h . 
these models assume a fully sampled stellar IMF, which 
for smaller mass clusters could lead to the inclusion of 
unphysical fractions of stars. With stellar clusters of 
high masses (> 1 x 10'^ Mp^), we expect the IMF to b e 
fully populated (jElmegreenl 120061: Ida Silva et all 120121 ) . 
so this is not a problem, but, as the cluster masses de- 
crease, massive stars are less likely to be formed and 
massive stellar populations are not fully represented. In 
order to properly measure parameters for low mass clus- 
ters (^^500 - 5000 M0), it is then important to turn to 
stochastic modeling. For this paper we will mainly fo- 
cus o n the stochastic models of SLUG (|da Silva et al.l 
120121 Stochastically Lighting Up Galaxies) which per- 
forms the synthesis of composite populations using in- 
dividual stellar clusters which are stochastically popu- 
lated with stars, using the IMF as a probability distri- 
bu tion function. Othe r stoch astic models are presented 
in IPopescu fc HansonI (I2010D, but are not used in this 
study. According to rCerviho fc Luridianal ([2006, 2004) 
clusters with masses below 10'^ M0 and ages less than 10^ 
years may be susceptible to color biases from determin- 
istic stellar synthesis models. NGC 4214 has numerous 
small clusters making it pertinent that stochastic mod- 
els be used. As a check and comparison, we also employ 
the deterministic SB99 models using the same input pa- 
rameters and use outputs which c ontain both s t ellar a nd 
nebular emission as cautioned bv lReines et al.l H2010D . 

Both sets of mo dels use a Kroupa IMF between 
0.08-120 Mq (jKroupa 2001), Padova AGB tracks with 
z=0.004, and assume that the clusters form in a single 
instantaneous burst. For the truncated SLUG models, 
where the maximum mass was only allowed to be 30 
M0, a Salpeter IMF was used. The SB99 models contain 
ages between 1-200 Myr in time steps of every 1 Myr for 
models between 1-20 Myr, and every 25 Myr for models 
between 25-200 Myr (to easily distinguish old clusters.) 
New SB99 models were generated between 2-8 Myr with 
time steps of 0.2 Myr for more accurate comparison with 
SLUG models. The SLUG models include about 40000 
cluster templates with ages from .01 Myr to 20 Myr with 
a cluster mass of 1 x 10^ Mq (Figure [2]). Tests have 
been run on a subset of clusters using 5 x 10^, 1 x lO'^, 
and 3 x 10^ M© SLUG models and we have found that 
the differences in inferred ages and masses of the same 
cluster among the various models is small, and is already 
encompassed by the uncertainties generated within one 
1 X 10^ M0 model, which is described below. Therefore, 
to simplify the analysis and most accurately reflect the 
masses of the clusters in NGC4214, we use only 1 x 10"^ 
M0 SLUG models for our cluster sample. SLU G models 
do not allow for binarity, but a recent study bv lEldridgi 
(|2012D indicates that at masses > 10^ M©, the scatter 
in Lfl-c/Mci is the same between models that use single 
stars and those that introduce binaries, so this should 
not introduce additional uncertainties. 

3.2. Age and Mass Determination 

To estimate the age, mass, and extinction of each clus- 
ter we employ a reduced fitting technique between 



both the SLUG and SB99 models and the cluster pho- 
tometry. To do this, we have used the Yafit (Yet Another 
Fitting ToolE3) curve fitting tool, which was created to 
fit photometry with model SEDs. It provides both a 
graphical and a numerical output containing the reduced 
value and the scaling factor between model and obser- 
vation. The ensemble of photometry for each cluster was 
compared to both SLUG and SB99 models spanning the 
complete reddening range between 0< E{B — V) <0.40. 
Note, we do not use the Ha filter as part of the fit, as we 
do not want to bias our sample based on the presence or 
absence of Ha emission. 

After all of the observations were compared against 
each model, the fits for each cluster were then organized 
into increasing reduced values. It is important here 
to point out that there is no single solution for the age 
and extinction of the cluster, but instead there is a range 
of best fits which could produce the model fit. In Figure 
[3] we show both the SED fits and the range of ages and 
extinctions spanned by the model fits to the youngest (~ 
2 Myr) and oldest (~ 7.5 Myr) cluster in our sample. 
By using the large range of ages and extinction consis- 
tent with the model fits, the actual mass distribution 
has been spread over a large range of values and only 
produces a peak at the most probable value. Therefore, 
we allowed the values t o rang e between and 1, as 
was done in iPasquali et al.l (|200>f ). and include all ages, 
extinctions, and therefore corresponding masses within 
that range. Masses were determined using the scaling 
factor output by Yafit and attributing error bars consis- 
tent wit h the error bars in the ages and extinction values. 
As with iCalzetti et al.l (|2010( ) , we also find that chang- 
ing the upper mass limit of the IMF in both the SLUG 
and SB99 models from 120 to 30 Mq increases the mass 
estimates by roughly 2.5. To double check our initial re- 
sults, comparisons with previous studies which have de- 
termined ages and masses of some of the clusters in NGC 
4214 were conducted and found to be comparable. For 
example, the young massive cluster located at the center 
of NGC4214 (R.A.=12''15"39^44, Dec = 36°19'34".94), 
noted here as Cluster 1, has age and mass estimates of 
4-5 Myr and 2.7±0.4 x 10"* Mq ()MacKentv et al.l 120001: 
iLeitherer et al.l [r996( ) using a Salpeter IMF between 1- 
100 Mq. Our best estimate, from the best fit from both 
SLUG and SB99, is 4.2 and 4.8 Myr, respectively, with 
a corresponding best fit mass of 9.3±4 and 9.7±3 x 
10^ Mq, which is consistent with these previous studies 
within 2(7. Unfortunately the cluster has blown an exten- 
sive asymmetric bubble, clearing the surrounding region 
of much of its hydrogen gas and allowing more than ha lf 
of the ionizing photons to escape (jMacKentv et al.l l200Cf) , 
making accurate measurements of LiHa extremely diffi- 
cult. For completeness purpose s, we will use the Ha flux 
for Cluster 1 (LAs) quoted in iMacKentv et al.l (2000,) . 
which gives an Ha luminosity of 8.4 x 10^^ erg s~^. 

By < 8 Myr, the compact HII region surrounding stel- 
lar clusters has expan ded into a shell struct ure which 
disperses into the ISM (iWhitmore et al.ll2011| ) and mas- 
sive stars capable of producing ionizing photons (> 15-20 
Mq) have disappeared. For example, SB99 models using 
the parameters listed above, show that the Ha luminos- 

http://www.star.bris.ac.uk/ mbt/yafit/ 
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Fig. 3 — Top: SED model fits for a ~2 Myr cluster (left, cluster 402 in Table 1) and a --7.5 Myr cluster (right, cluster 
188 in Table 1). Each fit uses an E{B — V) of 0.06, and 0.40, respectively. SLUG models are indicated by a solid line, 
and SB99 models by a dashed line. The large uncertainties in the photometry can allow for a large range of fits with an 
acceptable value. Bottom: Contours for the reduced values for various age and attenuation values for a ~2 Myr 
cluster (left) and a '^7.5 Myr cluster (right). 

ity at 8 Myr is only 2.5% of Lho at 2 Myr. Therefore 
it is necessary to exclude all clusters > 8 Myr in order 
to keep objects in which the HII region is still density 
bound and the massive stars are still retained. This en- 
sures that uncertainties in the ionizing photon rate are 
reduced and that constraints can still be made on the 
upper end of the IMF. As will be discussed below and 
is shown in Figure |4l including ages greater than 6 Myr 
may already be too old for this type of analysis. We must 
also be careful of confusion, for example, including mul- 
tiple objects which may share the same HII region. For 
this reason, we have excluded those clusters which may 
share in ionizing common gas, removing clusters which 
are in excessively crowded regions which would hinder 



the correct measurement of I^Ha ■ This has left us with 
a total of 89 clusters between the ages of 2-8 Myr, 52 of 
which have masses that are > 500 Mq , and either have a 
measured Ha luminosity (38) or have an Ha luminosity 
that is non-detectable down to the 3(t limit of 1.6 x 10^^ 
erg s~^. Ten of these objects (Table 1, footnote c) have 
PSFs that are consistent with a single massive star, yet 
their SEDs require the presence of multiple stars to be 
fully accounted for in flux. These low-multiplicity clus- 
ters tend to be among our lowest mass systems (Table 
1), and may be the NGC 4214 equivalent of the Trapez- 
ium cluster in Orion, which has a 1-1.5 pc size and a 
handful of 15-30 M0 stars. Therefore, they have still 
been included in the sample and are indicated in Table 
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1 with the rest of the clusters. Table 1 also includes the 
corresponding ages, extinctions, and masses from both 
SLUG and SB99 for aU 52 clusters; Figure [T] shows their 
placement in the galaxy. We hav e also included the clus - 
ter naming nomenclature used in M acKentv et al.l ()2000f ) 
where applicable for ease of comparison. 

4. DISCUSSION 

If the IMF is populated purely stochastically, then we 
can expect that 100 10^ Mq clusters would contain the 
same numbers and masses of stars as one 10^ Mq clus- 
ter, an d that both would represent a fully sampled IMF 
(|Elmegrc en 2001. 2006). On the other hand, if the most 
massive star in a cluster is limited by the mass of the 
parent cluster we might not expect any high mass stars 
to be present in 1000 Mq clusters. For e xample, in the 
variable upper mass limit formulation of IWeidner et afl 
dJOlO), they propose a M(max)* - Mc; relation in which 
no stars more massive than 35 Mq would be present in 
a lO'^ Mq cluster. The summation of the total ionizing 
flux from these small clusters divided by the total cluster 
mass should then be much lower than the ionizing flux 
from the single large cluster divided by its mass, and as 
cluster mass decreases there is a deviation from the ra- 
tio of ionizing photons to mass expected by a universal 
IMF (Figure [SJ dashed-dotted line) . In a universal IMF 
scenario, then, the summation of the total ionizing flux 
from the small clusters divided by the mass will be con- 
sistent with that of a single large cluster. Even though 
most low mass (< 500 Mq) clusters will produce low Ha 
luminosities in the universal scenario, there will be some 
that do produce a larg e ionizing continuum fro m the odd 
star well over 20 Mq. IVillaverde et"all (IMol) estimates 
that only 20% of 100 Mq clusters will have stars large 
enough to create an H II region, therefore given a large 
enough sample the effects will average out. By summing 
the Lho and masses of all of the small clusters into one 
data point, not only are the observational uncertainties 
reduced but the stochastic effects are minimized. 

We have used three mass bins (see Figure [S]), each with 
a mean mass of 7.5 x 10^ Mq, 2.2 x 10^ Mq, and 4 x 
10^ Mq. The error bars in each bin are obtained by 
adding in quadrature the individual mass and luminos- 
ity errors of each cluster fit. The largest mass bin has 
only 5 members in total, including the super star cluster 
at the center of the galaxy (Cluster 1), for which the Ha 
luminosity i s highl y uncertain and we are using a value 
obtained in iMacKentv et al.l (|200g) . These extremely 
large clusters with negligible Ha emission are included 
to illustrate the possibility that we may have feedback 
occurring in our larger clusters which is dispersing the 
gas more efficiently. Indeed we note that the most mas- 
sive star clusters in NGC 4214 are surrounded by ionized 
gas shells (see below). We have also included for illus- 
tration those clusters that have masses greater than 500 
Mq, but only upper limit measurements for Lhu (1-6 x 
10^^ erg s~^, open squares). Finally, we have included 
the full sample of objects with masses greater than 500 
Mq with those objects which may be single stars or low 
metallically clusters removed. In the lower mass bins, 
the error is dominated by the range of masses. Also plot- 
ted in Figure [S] is the expected average Lna/^ci from a 
1/5 Zq SB99 metallicity model that is fully populated 
up to 120 Mq between 2-5 Myr (top dashed line), and 



2-8 Myr (bottom). The expected range for a M*- Md 
model where the most massive star in the cl uster is a 
function of cluster mass ([Weidner et al.|[2010[ ) is shown 
in the dashed-dotted line also averaged between 2-5 Myr 
(top) and 2-8 Myr (bottom). If the metalhcity of NGC 
4214 is 1/4 Zq we do expect the data to fall somewhat 
below the stellar synthesis models, since lower metallici- 
ties create higher Ha luminosities. 

When including upper limits in Ha the clusters in 
the highest mass bin have a value of hna/^ci consid- 
erably lower than that of lower mass clusters (Figure [5]) . 
We have investigated whether the clusters in the high- 
est mass bin may be experiencing feedback effects which 
have caused ex p ulsion of the gas from the HII region. 
iPellegrini et al.l ()2012l ) have found that the amount of 
ionizing photons lost from a cluster can be dependent 
on the HI density surrounding the cluster. With this 
in mind, we have used the H I maps of NGC 4214 pub- 
lished in I Walter et al.l ()2001|) to locate the clusters in the 
largest mass bin without Ha emission. We find that they 
are located in regions that have HI densities roughly 2/3 
that of the maximum density, which by itself would not 
indicate that the clusters should experience significant 
ionizing photon loss. When observing the clusters in the 
Ha image though, it is quite clear that some mechanism 
has blown much of the gas away from the clusters and 
they are surrounded by wind blown bubbles. This is not 
completely unexpected, as these clusters with extremely 
low L ho values are located in th e NGC 4214-1 region, 
which lMaiz-Apellaniz et ahl (|1998[ ) found there was a sig- 
nificant decoupling of the stellar clusters with the ionized 
gas. 

In Figure |4] we present histograms of the ages of the 
stellar clusters with detected Ha (blue) and non-detected 
Ha (red) in our sample, including one containing all clus- 
ters younger than 8 Myr with masses < 10^ Mq (left), 
and all clusters with masses between 500 - 10000 Mq 
(right). What we have found is that all clusters with ages 
that lie between 6-8 Myr have non-detections in Ha. All 
clusters younger than 4.5 Myr are detected in Ha, and 
2/3 of clusters between 4.5 and 6 Myr are detected in 
Ha. We must note that the uncertainties in each bin of 
Figure 4 are larger than the bin size, yet even when ac- 
counting for the uncertainties in the best fit age, the ages 
of clusters showing Ha remain below 6 Myr. The excep- 
tion here arc those few clusters with masses > 10^ Mq. 
These have not been included in Figure 2] because even 
at ages between 2-3 Myr, they have very little measured 
Ha flux, likely due to stronger feedback effects from the 
clusters themselves, or the more crowded environment in 
which they are located which makes it difficult to deter- 
mine which ionizin g photons come from which cluster. 
iRelafio et al.l ()2012D have also determined that the leak- 
age of ionizing photons is expected for those HII regions 
with ages greater than 4 Myr. While some of the de- 
tected clusters do venture into ages greater than 5 Myr, 
we do need to be aware of the fact that especially in clus- 
ters that may contain only one or two extremely massive 
stars they may not live long enough to produce ionizing 
photons out to 8 Myr. For example, the lifetime of a 35 
Mq star is roughly 5 Myr, while a 15 Mq star may live 
100 Myr and a 120 Mq star only 2 Myr. If there is only 
one massive star in these smaller clusters, the L^f^ may 
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Fig. 4. — Histograms of relative frequencies of best fit cluster ages for all young clusters (left) and clusters with masses > 500 Mq (right). 
The blue histogram indicates those clusters with measured Ha emission, while the red indicates those with only upper limits. From both 
plots it is apparent that clusters < 4.5 Myr will still retain Ha emission, while those > 6 Myr will not. Cluster with ages between 4.5-6 
Myr have a lower probability of producing Ha emission. These plots exclude those clusters with masses > 10* Mq as they may not be 
fully explored by the SLUG models and are far more likely to have expelled their surrounding hydrogen gas at a much younger age. 

be more sensitive to the age. 

For each model, we show a histogram (Figure 6) of the 
most massive star in the model which best describes the 
photometry in each of the 47 clusters located in the two 
least mass bins. Included in this plot are also the error 
bars associated with the ensemble of fits within Ax^ of 
10% of the lowest value, and model lines for a fully sam- 
pled IMF (red) and a model truncated at 30 Mq (blue). 
What we find is that 25% of these model clusters have 
a maximum mass between 20-30 Mq, 50% have a mass 
< 40 Mq, and the remaining 50% contain most mas- 
sive stars greater than 40 Mq. We must point out that 
because clusters less than 500 Mq were excluded from 
the study it is possible that we may be underestimating 
the m ost massive st a rs wit h masses < 20 Mq. Accord- 
ing to iKroupa et al.l (|2011[ ). specifically Figures 2 and 5, 
the expected maximum stellar mass of a c luster between 
500-30 00 Mq is between 20-40 Mq if the IWeidner et all 
([2010') hypothesis holds. Some of our clusters have a 
maximum stellar mass at this cluster mass in Figure [6l 
but the existence of stellar masses greater than 40 Mq 
for the same range of cluster masses indicates that there 
is no maximum mass for the stars in these clusters, other 
than the usual upper limit found in much more massive 
clusters. Our results in Figure [6l a re consistent with th e 
random sampling from Figure 2 of IKroupa et al.l (|2011l ). 
where the maximum stellar mass can range between 15- 
120 Mq for a 500-3000 Mq cluster. We must point out 
that these models were populated specifically for 1 x 10^ 
Mq clusters, and the clusters themselves range between 
500-9000 Mq, but that our treatment should give a rea- 
sonable representation of the actual spread of massive 
stars. As mentioned above, tests run on 500 Mq, 1000 
Mq, and 3000 Mq SLUG models resulted in mass and 
age discrepancies within the uncertainties of the < 1 
resu lts. 

Larsenl (|2002[ ) has suggested that lower SFR corre- 
lates with a lower average cluster mass; therefore most 
galaxies with a low SFR should be deficient in high 
mass clusters in a purely stochastic case. For exam- 



ple, 100 galaxies with a SFR of 0.1 Mq yr ^ will have 
the same distribution of cluster masses as 1 galaxy 



with a SFR of 10 Mq yr^^, yet most of the low SFR 
galaxies will be lacking massive clusters and a few will 
have an over abundance. Other studies have indi- 
cated that there is a systemat ic trend for small clusters 
to form only low nia,ss stars (IKroupa h Weidneii [20031 : 
IWeidner et all [20T0I : iPflamm-Altenburg et al.l I2009L for 
examp le) , effec tively steepening the the slope of the 
uIMF (jMeurer et al. .2009i) . When interpolated, across a 
galaxy, the integrated gal actic IMF (jKroupa fc Weidneii 
l2003HWeidner et al.ll20ia IGIMF) would ultimately lead 
to low SFR galaxies such as dwarfs to be deficient in both 
large stellar clusters and massive stars. From what we 
have presented above, we do not find that this is the case 
in NGC 4214. The SLUG models suggest that roughly 
50% of clusters with masses around 1 x 10^ Mq have at 
least one star more massive than 40 Mq (Figure [6]) . We 
have also performed a Komogorov-Smirnov (K-S) test be- 
tween the best-fit Ljjq l^ci data and both SLUG models 
sampled up to 30 Mq and 120 Mq (blue and red lines 
in Figure [71 respectively). We find that there is about a 
65% chance the data were drawn from the parent model 
of the fully populated Kroupa IMF up to 120 Mq, but 
only about a 1% chance they come from the models trun- 
cated at 30 Mq . By adding in error bars from the ensem- 
ble of fits within Ax^ (right), these numbers become 75% 
and 1% respectively . When we compare the Ljjq /Mci 
values of the clusters which were age-dated using a SB99 
model where the most massive star is only allowed to be 
30 Mq against the fully populated and truncated SLUG 
models (see Figure [71 bottom) the K-S test only shows 
agreement of 7% and 4% respectively that the data come 
from those respective populations, even when uncertain- 
ties are included. Thus, even allowing our clusters to be 
modeled as drawn from a parent population that does 
not have stars more massive than 30 Mq (which affect 
the determination of the cluster masses) the l^ci 
distribution does not change in such a way to agree with 
a truncated IMF. 

5. CONCLUSIONS 

We have attempted to constrain the upper end of the 
IMF of the nearby irregular star bursting dwarf galaxy 
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would result in the maximum stellar ma ss in a cluster 
of 10'^ M0 being no greater than 35 M0 (iWeidner et alj 
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Fig. 5. — Top: Location of Lhq, /M^i for mass bins in NGC 4214. 
Filled blue squares only include clusters with Ho emission (all of 
which are less than 6 Myr), while open blue squares include all clus- 
ters, and blue stars include all clusters without those objects which 
could possibly be single stars or low-multiplicty clusters. Dotted 
line is the expected L^f^ /M^j for a universal I MF while the dash- 
dotted line is for a variable upper mass limit IIKroupa fc Weid neil 
[2003: Weidncr et al. 2010) where the most massive star in a cluster 
is determined by cluster mass. The top lines show averaged models 
between 2-5 Myr, while the bottom is between 2-8 Myr. Bottom: 
Same as top, but also containing cluster measure ments of NGC 
4214 from SB99 (cyan circles), from M51 ( Calzctt TeFall I20T0I . 
red triangl es), R136a l[Pel lc grini et al. (2012), green diamond) and 
NGC 330 llPellegrini et aEI (goQ), purple diamond). The location 
of the largest mass blue and cyan symbols is highly impacted by 
the uncertainty of the Ha luminosity of Cluster 1. 



NGC 4214 using the methods of ICalzetti et all (|2010D . 
which uses the ratio of the luminosity of the ionizing 
photons normalized to the mass of the cluster as a proxy 
for probing the presence of massive stars. With a final 
sample of 52 young clusters with masses > 500 Mq, we 
have determined that even at masses ~ 10^ Mq, there 
does not seem to be a deviation from the expected ion- 
izing flux of an universal IMF up to 120 Mq. Clusters 
with a mean mass down to 700 Mq have a Lho /Mc; 
ratio that lies along that predicted by a universal IMF. 

We have also determined that a truncated IMF (one 
in which the maximum M* is a function of Mcj), which 



0.3500 r 




0.0500 - 
0.0000 U 

' I ' ' 

40 60 80 100 
Mass of Largest Star in Cluster 

Fig. 6. — Histogram of the most massive star in each model 
considered the best fit along with fits within Ax^ of 10 % of the 
lowest value from SLUG plotted against a fully sampled IMF (red) 
and a model truncated at 30 Mq (blue). Masses range from 10-119 
Mq , and while 50% of the most massive stars are < 40 Mq , there 
are still the existence of stars between 40-120 Mq in clusters with 
masses of 10^ Mq which would be contrary to a variab le upper 
mass limit I IKroupa fc Weidneril2003 : IWeidner eraLll2010l ) theory. 



I2010f) . does not sufficiently explain the young clusters 
seen in NGC 4214. Models used to age-date the clusters 
indicate that up to 50% of the clusters contain a mas- 
sive star greater than 40 Mq , while Komogorov-Smirnov 
tests indicate that there is only a 1% chance that the 
L_f/a /Mc; values from clusters in NGC 4214 come from 
a parent model with the maximum mass truncated at 
30 Mq and a 75% chance that they a rise from a fully 
popul ated IMF up to 120 Mq. As with iFumagalli et alJ 
(|2011f ). who investigated the integrated properties of in- 
dividual galaxies, we find that the summation of indi- 
vidual young clusters in NGC 4214 is more consistent 
with a universal IMF without a truncation of massive 
stars. Our test shows not o nly that the M, - Md rela- 
tion of IWeidner eTall (poTol ) is in disagreement with the 
observations, but also that the data are consistent with 
a stochastically-sampled I MF with our upper mass limit 
consistent with a standard iKroupal (|2001i) IMF. 
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Fig. 7.— Histo grams of \^Ha /^cZ from SLUG models for a fully sampled IMF (red) and a model truncated at 30 IVIq (blue) for Lj/ck 
/Mcj plotted against clusters from NGC 4214 whose masses were determined from SLUG models with a maximum stellar mass of 120 Mq 
(top) and SB99 models with a maximum stellar mass of 30 Mq (bottom). The left panels only include the single best-fit mass, while the 
right panels include error bars for all masses with a < 1- 



APPENDIX 
FILTER CONVOLUTION 

SLUG models, as of the submission of this paper, are not delivered in HST filter pass bands. Luminosities were 
delivered in Johnson- Cousins UBVI and Galex NUV, and, because of this, transformation coefficients were derived 
to approximate WFC3 filters. To derive these coefficients, each 1 Myr time step of the SB99 model was convolved 
with the appropriate WFC3 filter passband (discussed below) to get the fluxes in each filter at each age. The SLUG 
models were then divided into 1 Myr age bins (for example, the 3 Myr age bin consisted of all models between the 
ages of 2.5 and 3.4 Myr), and the fltxx in each filter was averaged into a mean flux. Therefore for 1 Myr age bins, 
there exists an average luminosity for each filter to directly compare with SB99 models. The ratio between the SB99 
and SLUG models was then taken, resulting in coefficients to convert SLUG filters to WFC3 filters. These values were 
then applied back onto each of the individual ^ 40000 SLUG models, and the effective wavelengths were adjusted to 
refiect those of WFC3 filters. 

For the SB99 models, the filter convolution process was more straightforward. From the output SED each f Myr age 
bin was extracted and then using the IRAF function sinterp were all interpolated into the same wavelength range of 
500-10000 A. Each WFC3/UVIS filter transmission curve was also extracted into the identical wavelength range, and 
then integrated over those wavelengths to determine the total transmission for each filter at the effective wavelength 
(/ T\d\). The SB99 SEDs were then multiplied by the transmission curves and a numerical integrated value was 

obtained (J F\T\dX). Next, < F\ > = ^ J^^^'J^ is then calculated for each age bin and each filter. This transforms 

the continuous SB99 SED into 5 distinct photometry points which can be easily compared with the photometry of the 
NGC 4214 clusters (Figure 3). 

Once the models have been convolved with the filters, they are then corrected for possible host galaxy extinction . 
Due to the low metallically of this galaxy, we have chosen to use an SMC extinction curve from iFitzpatrickl (|19990 
which is more in line with the metallically of NGC 4214. We have limited our reddening to lie between 0.0 < E[B — V) 
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< 0.40 as previous studies on this galaxy seem to indicate low values of extinction (|Ubeda et al.ll200"7f) . The population 
synthesis SEDs are first convolved with the extinction curve at selected values of the color excess E(B- V), spanning the 
full range of E(B-V) used for our analysis (0.0-0.4), order to determine the effective wavelengths Xext,eff to be used 
when applying extinction to broad band filter photometry. In general Xext,eff is different from the effective wavelength 
of a filter, owing to the non-symmetric transmission curves of most filters. We should note that Xext.eff depends also 
on the range of E(B-V) considered, and should be re-calculated when using color excesses outside our range. We elect 
to apply extinction corrections after filter convolution because tests comparing the results of this procedure against 
the other procedure of applying the extinction before filter convolution gives differences of only 5%. Our choice, 
however, provides the flexibility to change E{B — V) at will. Each filter is then multiplied by ^'^■^'^^e(b~v)^SMc ^ 
where E{B — V) is in increments of 0.02, and the appropriate SMC extinction value is determined for the wavelength 
of each filter. 
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ID a (J2000) 



5 (J2000) 

36° + 



Age 
(Myr) 



SLUG 
Mass 
(10' Mq) 



E(B - V) 


Age 


Mass 


E{B - V) 




Other 




(Myr) 


(10^ Mm) 




erg s"'^ 


Names^ 


0.24 


4.0 


0.21 


0.16 


1.6E+35 




0.16 


5.0 


0.12 


0.2 


1.6E+35 




0.32 


4.0 


0.50 


0.4 


1.6E+35 




0.22 


5.0 


0.33 


0.22 


1.6E+35 




0.22 


5.4 


0.20 


0.12 


1.6E+35 




0.3 


5.6 


0.20 


0.26 


1.6E+35 




0.4 


7.6 


0.40 


0.34 


1.6E+35 




0.32 


4.0 


1.00 


0.26 


1.6E+35 




0.06 


5.8 


0.49 


0.02 


1.6E+35 




0.12 


5.4 


0.46 


0.1 


1.6E+35 




0.06 


4.8 


3.60 


0.04 


1.6E+35 


I-B2n 


0.08 


4.8 


4.58 


0.08 


1.6E+35 


I-Bs 


0.06 


5.2 


3.77 


0.06 


1.6E+35 




0.38 


4.0 


0.21 


0.3 


2.3E+36 




0.16 


3.8 


0.38 


0.16 


3.1E+37 


H-Cln 





4.6 


0.17 





3.7E+36 




0.2 


3.2 


1.01 


0.22 


8.8E+36 




0.28 


4.2 


0.39 


0.28 


4.0E+36 




0.24 


5.2 


0.19 


0.16 


2.3E+36 




0.14 


3.2 


0.39 


0.14 


5.1E+37 


I-A5n 


0.06 


3.6 


0.22 


0.08 


5.1E+36 


LCln 





5.2 


0.13 





7.0E+35 




0.22 


5.4 


0.10 


0.14 


l.lE+37 




0.04 


4.8 


0.13 


0.04 


9.7E+35 




0.18 


5.2 


0.86 


0.14 


7.3E+36 




0.22 


5.0 


0.62 


0.14 


3.9E+36 




0.2 


4.0 


0.35 


0.12 


2.8E+37 







5.4 


0.17 





3.6E+36 




0.24 


4.0 


0.38 


0.16 


l.lE+37 




0.04 


5.0 


0.34 


0.02 


9.5E+36 




0.18 


4.0 


0.95 


0.16 


7.5E+37 


I-A3n 


0.4 


4.0 


0.77 


0.38 


3.9E+37 


IXn 


0.4 


5.0 


0.58 


0.32 


2.2E+37 


II-A 


0.12 


5.0 


0.85 


0.1 


2.5E+37 


I-D2n 





4.6 


0.39 





2.1E+37 




0.14 


3.2 


0.78 


0.14 


3.9E+37 


Il-Cln 


0.12 


4.4 


3.39 


0.08 


1.6E+37 


I-Cln 


0.08 


5.0 


0.57 





2.5E+36 







6.2 


5.75 


0.08 


l.OE+38 


II-B 


0.4 


5.0 


0.90 


0.4 


8.6E+35 




0.38 


4.8 


0.43 


0.1 


2.2E+36 




0.02 


4.8 


0.57 


0.04 


1.6E+36 







4.8 


0.56 


0.22 


1.2E+37 




0.28 


5.8 


0.26 


0.2 


7.9E+35 




0.28 


4.6 


1.20 


0.26 


1.9E+37 




0.2 


5.6 


4.29 


0.16 


2.8E+36 


I-B2n 





4.6 


1.24 


0.02 


1.9E+37 


II-C2n 


0.2 


4.8 


1.31 


0.2 


4.8E+37 


I-Gn 


0.18 


4.0 


3.22 


0.14 


5.1E+37 


II-B 


0.14 


4.0 


2.68 


0.12 


9.7E+37 


I-Aln 


0.12 


4.6 


7.03 
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^ Clusters above the line have only 3(t detections of Ha. Clusters 
1-4 likely contain Ha emission, but due to size and nearness of other 
clusters have blown cavities surrounding them. 

^ Other names come from MacKenty et al. 2000. 

" PSF consistent with a single star at the distance of NGC 4214 



(FWHM ~ 1.3 pc), and SED colors consistent with either a single 
O star with several low mass (B5 or later) stars or a few early B 
type stars. 

Luminosity derived from Table 3 in lMacKentv et al.l II2000I ). 



